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Abstract 

We will summarize results of calculations of the modes of oscillation 
trapped within the inner region of accretion disks by the strong-field grav- 
itational properties of a black hole (or a compact, weakly-magnetized neu- 
tron star). Their driving and damping will also be addressed. The focus 
will be on the most observable class: the analogue of internal gravity 
modes in stars. Their frequencies which corrrespond to the lowest mode 
numbers depend almost entirely upon only the mass and angular momen- 
tum of the black hole. Such a feature may have been detected in the 
X-ray power spectra of two galactic 'microquasars', allowing the angular 
momentum of the black hole to be determined in one case. 

PACS: 04.70.Bw; 97.60.Lf 
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1 Introduction 

In this review, dedicated to one of my oldest and best friends, Giora Shaviv, on 
the occasion of his 60th birthday, a potentially powerful probe of black holes 
will be described. But first we may ask to what extent their existence has been 
verified, keeping in mind our definition of a black hole: a region of spacetime 
described by the Kerr metric. There are several lines of evidence. 

A) The determinations that a certain amount of mass is contained within a 
certain radius (typically many orders of magnitude greater than the horizon size 
of a black hole of that mass) are based mainly upon (a) the time dependence 
of the Doppler velocity of a companion for binary systems and (b) the Doppler 
velocity and orbital radius of gas or stars for galactic nuclei. There is strong 
evidence for a mass greater than that of standard neutron stars in roughly six 
Galactic binaries (McClintock 1998), and of supermassive objects in roughly 
twice as many galactic nuclei. The rapid variability and high energies of much 
of the radiation provides evidence that the emission region is indeed compact 
(GM/Rc 2 > 0.1). 
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B) The energy dependence of emission line profiles, when interpreted as 
Doppler and gravitational shifts of excited line emission from accretion disks, 
provides values of its extent, in terms of R/M (e. g., Fabian et al. 1989; Tanaka 
et al. 1995). However, a dependence of the line emissivity on radius must be 
assumed. In one case, it has been claimed that a rapidly rotating black hole 
(a = cJ/GM 2 « 1 is required to explain inferred values of inner radius < 
2GM/c 2 (Iwasawa et al. 1996). 

C) Another approach involves fits to the energy dependence of the continuum 
(X-ray) spectrum (Zhang et al. 1997). A disk luminosity and maximum tem- 
perature is obtained for sources of known distance in which the non power-law 
component of the energy spectrum can be fit with a standard (diluted black- 
body) optically thick accretion disk model. This provides an effective area. 
After radiative transfer and general relativistic corrections, an inner radius of 
the disk is extracted. For luminosities sufficiently below the maximum (Edding- 
ton luminosity), is mainly a function of M and a. 

D) Recently, Narayan et al. (1997,1998) have claimed evidence for an event 
horizon, based upon the advectively dominated class of accretion disk models 
that have been developed (by Abramowicz, Lasota, Narayan, and others). These 
flows carry most of the energy generated into the black hole, rather than being 
radiated. Without a horizon, the energy would appear when the flow hit the 
surface of the central object (such as occurs with a neutron star). From spectral 
fits and the determination of the range of luminosities of systems whose accretion 
rate varies, they have identified various Galactic X-ray black hole candidates, 
as well as our and other galactic nuclei, as being in this state. An accretion 
rate less than about 8% of the maximum is required. However, there is still no 
criterion to determine whether a disk with such an accretion rate will instead 
choose to be in the standard (radiating) state. The low luminosity would then 
be due to a lower accretion rate. 

E) The frequencies of narrow features in the power spectra of various black 
hole candidates have been ascribed to persistent structures in the associated 
accretion disk. There have been various proposals for an identification with the 
orbital frequency of hot 'blobs', usually at the innermost radius of the disk. 
However, it is not clear (a) how such 'blobs' form, (b) how they survive as co- 
herent structures, and (c) why they are confined to a particular radius. Another 
identification is with inertial-acoustic traveling waves, which we shall comment 
upon later. Instead, we focus on the spectrum of normal modes of oscillation, 
which must exist at some level (determined by the driving and damping pro- 
cesses in the disk). In the same spirit with which helioseismology probes the 
interior of the Sun, this probe of the Kerr metric (and its accretion disk) has 
been dubbed (relativistic) diskoseismology. We now analyze this approach. [For 
an up-to-date survey of black hole accretion disk theory, see the monograph by 
Kato, Fukue, and Mineshiga (1998).] 

Since 1990, our group has been investigating consequences of the realization 
[initially by Kato and Fukue (1980)] that general relativity can trap normal 
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modes of oscillation near the inner edge of accretion disks around black holes. 
The strong gravitational fields that are required can also be produced by neutron 
stars that are sufficiently compact (requiring a soft equation of state) and weakly 
magnetized that there is a gap between the surface of the star and the innermost 
stable orbit of the accretion disk. Although we shall not explicitly consider such 
neutron stars here, the results obtained will also apply to them to first order 
in the dimensionless angular momentum parameter a = cJ/GM 2 , since their 
exterior metric is identical to that of a black hole to that order. It should be 
noted that a 0.5 for almost all models of rotating neutron stars (Friedman & 
Ipser 1992; Cook, Shapiro, & Tcukolsky 1994). 

These modes of oscillation provide a potentially powerful probe of both 
strong gravitational fields and the physics of accretion disks, since: 

• They do not exist in Newtonian gravity 

• Their frequencies depend upon the angular momentum as well as the mass 
of the black hole 

• The fractional frequency spread of each mode depends upon a, the elusive 
viscosity parameter of the accretion disk. 

2 Adiabatic Normal Mode Oscillations 

Following exploratory calculations by Kato & Fukuc (1980), Okasaki et al. 
(1987), and Kato (1989), we also employed a modified Newtonian potential to 
calculate the adiabatic eigenfunctions and eigenvalues of the lowest acoustic (p) 
modes (Nowak & Wagoner 1991) and internal gravity (g) modes (Nowak & Wag- 
oner 1992, 1993). This was extended to full general relativity by Perez (1993), 
who also included the corrugation (c) modes studied by Kato (1990, 1993) and 
Ipser (1994, 1996). The g-modes have been analyzed more extensively by Perez, 
Silbcrgleit, Wagoner, & Lehr (1997), and we shall also summarize their key re- 
sults below. 

We have applied the general relativistic perfect fluid perturbation formalism 
of Ipser and Lindblom (1992) to thin accretion disks in the Kerr metric. The 
radial and vertical components of the velocity of the fluid is neglected, which has 
little effect on the modes for these thin models. Neglecting the gravitational field 
of the disk (also a good approximation), the adiabatic oscillations of all physical 
quantities can be expressed in terms of a single scalar potential SV(r, z) cx 
Sp/p, governed by a second-order partial differential equation. The stationary 
(d/dt = 0) and axisymmetric (d/d<fi = 0) unperturbed accretion disk is taken 
to be described by the standard relativistic a-disk model (Novikov and Thorne 
1973, Page and Thorne 1974). We employ gravitational units (G = c = 1). 

All fluid perturbations are of the form f(r, z) exp[i(at + m<f>)]. With the 
fluid angular velocity, the comoving frequency is u>(r) = a + mf2(r). It is 
sufficient to consider eigenfrequencies a < and axial mode integers m > 0. 
Nonaxisymmetric modes (m ^ 0) should produce relatively little luminosity 
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modulation, unless the disk is viewed close to edge-on. In general, the vertical 
extent of the modes within the disk is restricted by the requirement that \u>\ be 
greater than the buoyancy frequency. Numerical results have been obtained for 
accretion disks which are barotropic (e. g., isentropic) on scales of order their 
thickness, in which case this restriction is not operative since the buoyancy 
frequency vanishes. 

The effective radial wavelengths A r are significantly smaller than r, allowing 
a WKB expansion of the problem and approximate separation of the governing 
equations, with SV — Vr(r)V v (r, 77). The vertical separability coordinate is 
77 = z/h(r), with h the effective half thickness of the disk. It and V v are slowly 
varying (compared to V r ) functions of r. 

The radial dependence of the fluid perturbations are governed by the equa- 
tion 



where Kerr metric component, U* is the fluid four-velocity component 

dt/dr, c s (r,z) is the speed of sound, and \&(r) is the slowly-varying separation 
function. The eigenfunction W(r) — (k 2 — oj 2 )^ 1 dV r / dr is proportional to the 
radial component of the fluid displacement. The key controlling ingredient in 
this problem is the relativistic behavior of the radial epicyclic frequency k{t), 
shown in Figure [l] with other important orbital frequencies. As usual, k is the 
frequency of radial perturbations of circular orbits of a free particle, while Q± 
is the same for vertical perturbations. 

From equation (|l]), we can see where the modes are trapped {W~ Y d 2 Wj dr 2 < 
0). Their classification is governed by the value of the eigenvalue as follows. 

• The g (inertial-gravity)-modes [defined by \& > (w/f2j_) 2 ] are trapped 
where uj 2 < k 2 . This occurs between some radii r_(u) and r + (c). The low- 
est modes (fewest number of radial and vertical nodes in their eigenfunctions) 
have approximately vertical displacements, and have eigenfrequencies \a\ = 2nf 
which are close to the maximum possible, a m — n(r m ) + md(r m ). The depen- 
dence of this frequency on black hole angular momentum is plotted in Figure ||, 
for m = 0. Plotted in Figure |§| is this radius r m to which r_ and r + converge 
as the frequency is raised to this maximum, as well as the effective radial width 
(r + — ?'_) of the eigenfunction. Although the approximate center of the eigen- 
function (located at r m ) is close to the temperature maximum of the disk, its 
relatively small width appears to limit the luminosity modulation to at most a 
few percent. 

The lowest radial (m = 0) g-modes have a frequency given by 



The properties of the accretion disk enter only through the small mode-dependent 
term e n j, which involves the disk thickness 2h(r m ) and the radial (n) and ver- 




(1) 
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tical (J) mode numbers, with S ~ 1. Typically h(r m )/r m ~ O.lL/LEdd for 
a radiation-pressure dominated optically thick disk region. F(o) is the upper 
function plotted in Figure || From Figure ||, we also see that the higher axial 
modes have a somewhat different dependence upon the angular momentum of 
the black hole, which in principle would allow its determination as well as that 
of the mass. 

• The p (inertial-acoustic)-modes [defined by if? < (ui/il±) 2 ] are trapped 
where u> 2 > k 2 . The radial (m = 0) p-modes that are trapped between the 
inner radius of the disk (r,) and r_ have very little radial extent, and thus 
will produce relatively little direct luminosity modulation, although they will 
modulate the accretion onto the black hole (or neutron star). The highest 
frequency modes that are trapped within r + and the outer radius of the disk 
(r ) will modulate a significant fraction of the disk where the luminosity per 
unit radius is highest. However, since their wavelength X(r) ~ 2ttc s /lo will 
be relatively short, their damping will be enhanced and the net luminosity 
modulation will be reduced. We have not yet investigated this outer branch of 
p-modes because of the uncertain physics and location of the outer radius. 

It should be noted that Honma, Matsumoto, & Kato (1992) and Milsom & 
Taam (1997) have claimed that inertial acoustic waves can be generated in disks 
with sufficient viscosity and fairly low mass accretion rates. These axisymmetric 
waves appear to propagate from the radius r m where k achieves its maximum, 
and have the corresponding frequency of the fundamental g-mode. However, 
the dependence of their wavelength A r = 2ir /k r on radius does not agree with 
the simple dispersion relation they assumed (lo 2 = k 2 + k 2 c 2 ). Possibly g-mode 
oscillations have coupled into these outgoing waves. 

• The c (corrugation)-modes [defined by "5? = (w/f2j_) 2 ] are typically nonra- 
dial (m = 1) vertically incompressible waves near the inner edge of the disk that 
slowly precess around the angular momentum of the black hole. The changing 
projected area of the mode could modulate the luminosity. The eigenfrequency 
of the fundamental mode and the corresponding trapping radius r c are given by 

2tt/ = |a| = n(r c ) - n±(r c ) , r c ^r t + CMiaT 1 ' 3 - 1) . (3) 

The mode extends from to r c . The constant C depends upon the speed of 
sound in the disk; a typical value is C = 0.17. Thus the mode has significant 
radial extent only for a <C 1, in which case the frequency becomes small and 
equal to the Lense-Thirring value: \a\ — > 2aM 2 /r 3 ,. The dependence of the 
frequency / on black hole angular momentum is also plotted in Figure |[ 

The frequencies of all of these modes are proportional to 1/M, but their de- 
pendence on the angular momentum of the black hole is quite different (c. f. Perez 
1993). It appears that the g-modes are the most robust, however, because their 
trapping does not involve the uncertain dissipative properties of the inner or 
outer disk radius. They should also be the most observable, because they usu- 
ally occupy the largest area of the disk, near the temperature maximum. 
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Also plotted in Figure || is the widely invoked frequency Vtmax/2-K of a 'blob' 
which orbits at the inner radius rj of the disk. Note that the ratio of this 
frequency to that of the fundamental g-mode remains fairly close to 3.08 (its 
value for a = 0) at larger values of a. 

3 Oscillation Amplitude and Frequency Width 

Nowak and Wagoner (1992) have estimated the damping (or growth) rates of p 
and g-modes due to isotropic and anisotropic turbulent viscosity and gravita- 
tional radiation reaction. Although the gravitational radiation instability will 
cause the modes to grow for M 3> 10 8 M e , we have found that isotropic viscosity 
will always produce damping. With j and n of order the number of vertical and 
radial nodes in any particular eigenfunction, its corresponding quality factor Q 
and frequency broadening is given by 

l/Q jn « (Af/f) jn ~ {j 2 + (h/r)n 2 ]a , (4) 

where the viscosity parameter is thought to lie in the range 10 -2 < a < 1. 
However, the effective width of the g-mode feature will also be determined by 
what range of mode numbers j and n are sufficiently excited. 

An estimation of the amplitude of the oscillations, and the resulting frac- 
tional modulation of the disk luminosity, was also obtained by Nowak and Wag- 
oner (1993). They found approximate equipartition between the energy of the 
lowest modes and the energy of the largest turbulent eddies. The turbulence 
can be driven by the Velikov-Chandrasekhar-Balbus-Hawley magnetic insta- 
bility [for relevant numerical results, see Stone et al. (1996)]. This produces 
maximum g-mode displacements of order the thickness of the disk, which should 
produce maximum luminosity modulations 8L U /L U ~ 10~ 2 , if the photon fre- 
quency v > h~ 1 kT(max). 

If turbulence drives these oscillations, it should also contribute to the fluc- 
tuation power in a broad range of frequencies. Nowak and Wagoner (1995) 
have computed the resulting power spectral density, and found that it could 
match that of the best-observed black hole candidate Nova Muscae (Miyamoto 
et al. 1993), but only at frequencies / > 10 Hz. The power is predicted to drop 
rapidly (oc /~ 5 ) above the frequency f2 ma:E /27r. 

Another potential excitation mechanism is 'negative radiation damping', 
produced by radiative losses in a radiation-pressure dominated region of the 
disk (Nowak et al. 1997). This can be understood from the dispersion relation 
for our modes, obtained from equation (Q) in the extreme WKB limit: 

uj 2 w n 2 ± c 2 s k 2 , (p-mode: +, g-mode: -) . (5) 

This expression reflects the fact that the pressure forces work against the effec- 
tive radial gravity for the g-modes, so the loss of pressure stengthens the mode. 
Again, the growth rate is proportional to a. 
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4 Observations 



These modes of oscillation should modulate the quasi-thermal emission from 
accretion disks in the ultraviolet for supermassive black holes and at soft x-ray 
energies for stellar mass black holes. We now discuss two black-hole candidate 
binary X-ray sources whose power spectra (obtained with the RXTE satellite) 
contained a feature which could be produced by a g-mode. The strongest ev- 
idence is the fact that (unlike lower frequency features that have been seen in 
various black hole candidates) the frequency did not change as the source lumi- 
nosity (and presumably mass accretion rate) did. This lack of dependence on 
the properties of the accretion disk is clearly predicted by equation (|J) for the 
lowest g-mode states. It is intriguing that these two sources are the strongest 
'microquasars' in our galaxy. 

Morgan, Remillard, and Greiner (1996) have detected a feature at / = 67 Hz 
with A/// k 0.05 in the microquasar GRS 1915+105. During the two periods 
that the feature was detectable, the source luminosity varied by a factor of ~ 2 
while the frequency changed by less than 3%. The amplitude of the peak varied 
from 1.5% at 3 keV to 6% at 15 keV. If this feature is produced by a g-mode 
oscillation, equation (|J) predicts a black hole mass of 10. 6Mq if it is nonrotating 
to 36.3M0 if it is maximally rotating. Other aspects of this identification are 
explored by Nowak et al. (1997). 

Remillard et al. (1997) have detected a feature at / = 300 Hz with A/// w 
0.4 in the microquasar GRO J1655-40. It was only detected during 7 observa- 
tions that showed the hardest X-ray spectra. Although the feature was clearly 
identified (with an amplitide of 0.8%) only by summing all observations, its 
frequency appeared to be stable. From extensive optical observations of the 
companion star, Orosz and Bailyn (1997) have determined the black hole can- 
didate mass to be 7.02 ± 0.22M Q , and the inclination of the binary to be 70°. 
This agrees with the estimated direction of the jet (if it is perpendicular to the 
orbital plane and the accretion disk) in both microquasars. 

In Figure |E|, we plot the relations between black hole angular momentum and 
mass if we identify the feature as being produced by the fundamental g-mode, 
c-mode, or 'blob' at the inner disk radius. Using a distance of 3.2 kpc and the 
observed energy spectrum of this source, Zhang, Cui, & Chen (1997) obtained 
an inner disk radius of approximately 20 km (by the method mentioned earlier). 
This yields the other relation between a and M shown. The value of the mass 
determined by Orosz & Bailyn (1997) is also indicated. 

Note that only by identifying the frequency with that of the g-mode is agree- 
ment with the spectroscopic determination attained, leading to a black hole an- 
gular momentum approximately 93% of its maximum possible value. If further 
tests confirm this identification, it will be the first accurate determination of 
the angular momentum of a black hole. 

One apparent problem is the fact that in both microquasars, the modulated 
flux has a far harder spectrum than that of the 1-2 keV blackbody emission 
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from the disk. However, the g-mode modulated luminosity will come from a 
radius significantly less than the average radius where the total luminosity is 
emitted. Consider then the effect of an atmosphere (or 'corona') of very hot 
electrons [temperature T a (r)] which scatters the photons emitted from radius 
r N(r) times. If the Compton parameter y((x T a N) is less than unity and 
increases as r decreases (as seems reasonable), then the Compton scattering can 
significantly harden the spectrum of the modulated photons relative to the total. 
Reasonable limits on the optical depth will prevent the photons from becoming 
demodulated by the electron scattering. 

5 Conclusions 

Returning to our initial question, we have seen that the evidence for compact, 
nonstcllar masses is strong, while the evidence for horizons is less strong but 
intriguing. The ultimate evidence of black holes (the Kerr metric) and measure- 
ment of their only other property (angular momentum) is still elusive. However, 
we have seen that the 'diskoseismic' probe, as well as spectroscopic (line and 
continuum) methods, have great promise. For instance, we note the fact that 
the coincidence between the g-mode and Stefan-Boltzmann determinations of 
the GRO J1655-40 black hole angular momentum would have disappeared if the 
metric was significantly different from Kerr. 

It is clear that continued X-ray timing observations at the times scales suf- 
ficiently short to see the high frequency modes are crucial to their identifica- 
tion. These are being carried out by RXTE (although power spectra to higher 
frequency are critical) and will soon be carried out by the USA satellite. In 
addition, we should not lose sight of the fact that supermassive black hole ac- 
cretion disks should exhibit similar signatures, merely scaled by the mass. For 
instance, the g-mode period for a 10 s solar mass slowly rotating black hole is 
1.6 days. Long-term monitoring with at least 1% accuracy photometry would 
be required. 

Most of the results reported here were obtained in collaboration with Dana 
Lchr, Michael Nowak, and Alexander Silbcrglcit. The research was supported 
in part by NASA through ATP grant NAG 5-3102 to R.V.W. and grant NAS8- 
39225 to Gravity Probe B. 
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Figure 1: The radial dependence of the square of the fundamental free-particle 
frequencies that govern the modes of the disk: Kcplcrian (O), and radial (k) 
and vertical epicyclic. Three values of the black hole angular momentum 

parameter a = cJ/GM 2 are chosen. 



Figure 2: The dependence of three characteristic disk frequencies on the dimcn- 
sionless angular momentum of the black hole. 



Figure 3: (a) The black hole angular momentum dependence of the radius r m to 
which r_ and r + converge as the g-mode eigenfrequency \a\ — ► a m . Also shown 
is the radius j-j of the inner edge of the disk, (b) The dependence of the fractional 
effective width of the lowest eigenfunction, Ar/r m = [r+(cr) — r-(a)]/r m , on 
the angular momentum of the black hole. The same values of m are chosen as 
in (a). The accretion disk model has a barotropic equation of state, and a speed 
of sound corresponding to a luminosity L = O.lL^dd from a radiation-pressure 
dominated disk. 



Figure 4: (a) The dependence of the maximum radial (m = 0) g-mode eigenfre- 
quency on the black hole angular momentum parameter a = cJ/GM 2 , relative 
to its value at a = 0. (b) The ratio of the maximum eigenfrequency of some 
higher m modes to that of the radial mode. 



Figure 5: Relations between black hole angular momentum and mass for the 
microquasar GRO J 1655-40. The three frequencies indicated are identified with 
the 300 Hz feature in its power spectrum. The value of inner radius was ob- 
tained by the spectroscopic (Stefan-Boltzmann) method of Zhang, Cui, & Chen 
(1997). The black hole mass, as determined by Orosz and Bailyn (1997), is also 
indicated. 
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